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STRONG RESPONSE OF THE VERY BROAD EMISSION 
LINE IN LUMINOUS RADIO-QUITE QUASAR PG 1416-129 
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ABSTRACT 

We report new spectroscopic observations performed in 2010 and 2011 for 
luminous radio-quite quasar PG 1416-129. Our new spectra with high quahty 
cover both H/3 and Ha regions, and show neghgible hne profile variation within 
a timescale of one year. The two spectra allow us to study the variability of 
the Balmer line profile by comparing the spectra with the previous ones taken 
at 10 and 20 years ago. By decomposing the broad Balmer emission lines into 
two Gaussian profiles, our spectral analysis suggests a strong response to the 
continuum level for the very broad component, and significant variations in both 
bulk blueshift velocity/FWHM and flux for the broad component. The new 
observations additionally indicate fiat Balmer decrements (i.e., too strong H/3 
emission) at the line wings, which is hard to be reproduced by recent optically 
thin models. With these observations we argue that a separate inner optically 
thin emission-line region might not be necessary in the object to reproduce the 
observed line profiles. 

Subject headings: galaxies: active — quasars: emission lines — quasars: individ- 
ual (PG 1416-129) 



1. INTRODUCTION 

Because of the limited spatial resolutions of the currently available instruments, the 
geometry and kinematics of broad-line region (BLR) in Active Galactic Nucleus (AGN) is 
still an important and unresolved problem. This poor understanding results in an uncertainty 
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for the virial coefficient / tliat depends not only on tlie kinematics, structure and orientation 
of the BLR, but also on the virial width estimator used. / is equal to 3 for a spherical isotropic 
velocity distribution, while / can be several times larger if the BLR can be described as a 
circular rotation disk. The poor determination of / finally causes an uncertainty in the 
estimation of the mass of central supermassive black hole (SMBH) of an AGN in the virial 
assumption (e.g., Onken et al. 2004; Kaspi et al. 2005; Marconi et al. 2008; Woo et al. 
2010). 

Recent spectral profile modelings indicate that the BLR is much more complex than 
a single virial component. In addition to the traditional broad component (hereafter BC, 
FWHM ~ 10^ km s~^), a very broad component (hereafter VBC, FWHM ~ 10"^ km s"^) is 
usually required to reproduce the observed line profiles of both and Ha in some AGNs 
(e.g., Sulentic et al. 2000; Netzer & Trakhtenbrot 2007; Hu et al. 2008; Marziani et al. 2009; 
Zhu et al. 2009; Zamfir et al. 2010). The physical interpretation of the two components 
is, however, still under debate at present. Previous studies have frequently suggested that 
the VBC likely arises in optically thin clouds that are closer to the central engine than the 
BC since 1980's (e.g.. Shields et al. 1995; Sulentic et al. 2000; Zamfir et al. 2010). By 
examining a large sample of quasars selected from SDSS, Hu et al. (2008) proposed that the 
BC (the intermediate-line region in Hu et al. 2008) can be explained by an inflow because 
it is systematically redshifted with respect to the quasar rest-frame. The comprehensive 
photoionization models calculated by Korista & Goad (2004) indicate that the VBC can arise 
from inner BLR (i.e., not a distinct emission-line region) where the Balmer line emission rate 
(and also responsivity) is low. 

Sulentic et al. (2000) detected a violent variabihty in both continuum and H^ emission 
fine in radio-quite luminous quasar PG 1416-129 (hereafter 2000 spectrum for short). By 
comparing it with the spectrum (hereafter 1990 spectrum for short) taken by Boroson & 
Green (1992, BG92), the 2000 spectrum shows a dramatic disappearance for the BC of the 
H/3 emission line. However, the VBC can be still clearly identified in the 2000 spectrum. 
Basing upon these observational facts, Sulentic et al. (2000) argued that a large fraction of 
the VBC is likely produced in optically thin clouds that are closer to the central source. 

In this letter, wc rc-obscrvcd the luminous quasar PG 1416-129 (R.A.=14'*:19'":03^.8, 
Dec=13 ° :10':44", J2000, z = 0.129280) at 10 years after the demise of its H/3 BC to examine 
the variability of the H/3 line profile. In addition to H/3, the Ha emission line is also covered 
in our new spectra. The observations and data reductions are described in §2. §3 presents 
the results and discussion. The A cold dark matter (ACDM) cosmology with parameters 
ho — 0.7, f^M = 0.3, and VIa = 0.7 (Spergel et al. 2007) is adopted throughout the letter. 
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2. SPECTROSCOPIC OBSERVATIONS AND DATA REDUCTION 

Two new optical spectra were obtained by the National Astronomical Observatories, 
Chinese Academy of Sciences (NAOC) 2.16m telescope in Xinglong Observatory at January 
25, 2010, and March 29, 2011. The spectroscopic observations were carried out with the OMR 
spectrograph equipped with a back-illuminated SPEC 1340 x 400 CCD as the detector. The 
used grating is 600g mm~^, and the slit oriented in the south-north direction corresponds a 
width of 2."0. This setup finally translates into a spectral resolution of ~ 4.5A as measured 
from the sky emission lines and comparison arcs. This spectral resolution is comparable 
with the previous observations in BG92. The blazed wavelength was fixed at 6500A in 
both observations in order to allow the observed spectra cover the wavelength range from 
4600A to 6900A in the observer frame (i.e., cover both H/3 and Ha emission regions). The 
observations were carried out as close to median as possible. In each of the observations, the 
object was observed successively twice with the exposure time of each frame of 3600s. The 
two frames were combined prior to the extraction to enhance the signal-to-noise ratio and to 
eliminate the contamination of cosmic-rays easily. The wavelength calibration was carried 
out by the helium-neon-argon comparison arc taken between the two successive frames, i.e., 
at the position being nearly identical to that of the object. The observed flux was calibrated 
by the Kitt Peak National Observatory (KPNO) standard stars HD 117880 and Feige56 in 
the 2010 and 2011 observations, respectively (Massey et al. 1988). 

The 2-dimensional spectra were reduced by the standard procedures through the IRAF 
packag^, including bias subtraction, flat-field correction and cosmical-rays removal before 
the extraction of the 1-dimensional spectra. Each extracted spectrum was then calibrated 
in wavelength and flux by the corresponding comparison arc and standard. In the observer 
frame, the redshifted Ha emission profile is strongly distorted by the A-band telluric ab- 
sorption feature at AA7600-7630 due to O2 molecules. In order to properly model the Ha 
profile, the two telluric features around A6800 and A7600 were removed from each observed 
spectrum by the corresponding standard. The Galactic extinction was corrected for both 
spectra by the color excess E{B — V) taken from the NASA/IAPC Extragalactic Database 
(NED), assuming the Ry = 3.1 extinction law of our galaxy (Cardelli et al. 1989). The 
spectra were then transformed to the rest frame, along with the correction of the relativity 
effect on the flux, according to the narrow peaks of H/3. Both reduced spectra (hereafter 
2010 and 2011 spectra) are displayed in Figure 1 by the bottom two curves. The top two 
curves show the 1990 and 2000 spectra for a comparison. 



^ IRAF is distributed by the National Optical Astronomical Observatories, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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For each of the four spectra, we modeled the continuum by a broken-powerlaw super- 
posed by a broadened Fe II template using a minimization. The template provided in 
BG92 is broadened to have a FWHM that is identical to that of total H/3 line. The mini- 
mization was carried out in the wavelength range free of strong emission lines (e.g., Ha, H/3, 
[O III]AA4959, 5007, He IIA4686) of AGNs. The continuum fittings are illustrated in Figure 
1 by the superposed red lines for the four spectra. 

The isolated emission-lines were then profiled by the SPECFIT task (Kriss 1994) in 
the IRAF package. We modeled each emission line by a sum of a set of several Gaussian 
components. As an illustration, the modelings in the 2011 spectrum are shown in the bottom- 
left panel and bottom-right panel in Figure 2 for the H/3 and Ha regions, respectively. The 
two middle panels show the same as the bottom ones but for the 2010 spectrum. In addition 
to the new spectra, the 1990 and 2000 spectra were also modeled by the same method in 
order to examine the variability of H^bc and H^vbc- The deblendings of the R/3 region are 
illustrated in the two top panels in Figure 2. 

3. RESULTS AND DISCUSSION 

The measured emission- line properties and continuum fluxes at 5 100 A are tabulated 
in Table 1. All the fluxes are reported after scaling the spectra to the total [O III]A5007 
flux of the 1990 spectrum, i.e., F([OIII]A5007) = (3.10 ± 0.23) x 10"^^ ergs cm-^ s'^ All 
the uncertainties presented in the table only include the errors resulted from the spectral 
modelings. 

3.1. Reliability of Spectral Analysis 

The red-shelf of the H^ hne proflle of the object is most unhkely contaminated by the 
broad He IAA4922, 5016 emission lines (Veron et al. 2002). At flrst, the broad He IA5876 
line emission is not detectable in both new spectra. Secondly, a corresponding red-shelf can 
be obviously identified in the Ha emission line because of the clear infiection of the profile. 
HavB is strongly redshifted with respect to the narrow core (e.g., Marziani et al. 2009) at 
a velocity shift At; ^ 2100 km s~^ in the 2010 spectrum, which is highly consistent with 
that of H/3vB, ^ 2400 km s~^. The consistence occurs not only in the velocity shift, but 
also in the hne width of ~ 10^ km s~^. Such consistence can be also identified in the 2011 
spectrum, although the Have bulk velocity is clearly less than that of H^vb (Note the large 
uncertainty for the H^ line). 
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We calculate the virial mass of the central SMBH (Mbh) from the four spectra taken 
in different epochs since we believe that the blackhole mass increases negligibly during the 
20 years. In each epoch, Mbh is estimated from Ha and H/3 emission lines by adopting the 
calibrations that are based on line width and line luminosity (e.g., Peterson et al. 2004; 
Kaspi et al. 2005; Green & Ho 2005; Bentz et al. 2006; Vestergaard & Peterson 2006). 
The luminosities are obtained after the correction of the local extinction. The extinction is 
inferred from the measured narrow-line flux ratio Hq;/H/3 by assuming the standard Case B 
recombination and a Galactic extinction curve with Ry = 3.1. The line width measured from 
the total broad-line proflle is usually adopted as a good virial broadening parameter because 
of the produced close virial relationship (e.g., Peterson & Wandel 2000; Onken & Peterson 
2002; KoUatschnny 2003;). By re-analyzing the reverberation- mapping data, Peterson et al. 
(2004) claimed that the best line-width measurement is the line dispersion calculated from 
the variable part of the spectrum. In order to measure the FWHMs of the entire broad-line 
profiles, the modeled narrow emission-lines are subtracted from each of the observed spectra 
before the measurements. In addition to the entire line profiles, Mbh is also estimated by 
using the FWHMs from the BGs. Note that the Mbh estimation is seriously complicated 
by the significant variation in the fine width and velocity shift between 2000 and 2010/11 
when the BGs are considered only. The VBGs are usually excluded from the Mbh estimation 
because of the large relative velocity with respect to the systematic velocity. 

Mbh is calculated from the fine widths and luminosities. The labels V and G are used 
to denote the calibrations in Green & Ho (2005) and Vestergaard & Peterson (2006), re- 
spectively. The calibrations finally yield consistent values of Mbh- For H/3 line, the average 
estimates are {M^^, Mb]^)bc+vbc =(2.5±1.3 x 10^ 1.8±0.8 x 1O^)M0, and {M^i, M^^)bc 
= (1.2 ±0.1 X 10^ 9.2 ±0.9 X 1O^)M0, where the formal errors are estimated from the obser- 
vations at the different epochs. For Ha line, we have average estimates of (M^^)bc+vbc= 
2.8toi X and {M^{^)Bc=8.2toi x lO'^ Mq from the two new observations. Note 

that these values are consistent with each other because an uncertainty of a factor of 4-5 is 
produced by the current available scahng relationships (Peterson 2008). 

3.2. Variations of BC and VBC 

One can see clearly that the object shows negligible line profile variation within the time 
scale of one year according to our new observations. The light-travel time cross a spherically 

1/3 —1/3 2/3 

symmetric BLR can be evaluated as tlt ~ 500Ljj^ 42/-7 days, where Lh/3,42 is the H/5 
fine luminosity in unit of 10^^ ergs s~^, /_7 is the fill factor of the emission gas in unit of 
10~^, and nil is the number density of the gas in unit of lO^^cm"^. The local-extinction 
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corrected luminosities are Lh/3 ~ 1.3 x 10^^ ergs s~^ for both BC and VBC components. 
The luminosities yield a light-travel time tlt ~ 500 day for both emission-line regions, which 
agrees with the observed fact that the line profile is non- variable within one year. 

The nature of the region emitting H^vb is still under debate. Many authors believe 
that the H/3vb components detected in many AGNs come from a separate inner optically 
thin emitting region according to the analysis of line ratios and the line-continuum variation 
studies (e.g., Fcrland et al. 1990; Marziani & Sulentic 1993; Corbin 1995, 1997; Corbin & 
Smith 2000; Sulentic et al. 2000; Zhu et al. 2009). The spectroscopic monitors previously 
claimed that the variations in line cores arc stronger than in wings for both high-ionization 
and low-ionization emission lines (e.g., Peterson et al. 1993; Kassebaum et al. 1997; Corbin 
& Smith 2000), although the opposite cases are found in Lya and C IVA1549 lines for high-z 
quasars (e.g., O'Brien et al. 1989). In an optically thin cloud, the emissivity of recombination 
lines is mainly determined by the volume and the density of the cloud rather than the 
incident ionizing continuum, which results in a weak response to the continuum level for the 
line intensity. By contrast, modern photoionization model of AGN's BLR suggests that the 
H/3vB component can be attributed to the inner part of the traditional BLR, because the 
calculated responsivity of Balmer emission of optically thick clouds increases steeply with 
the distance from the central engine (Korista & Goad 2004). Snedden & Gaskell (2007) 
compared their photoionization models with the broad emission-hne profiles of both HST 
and optical spectra. Their analysis points out that the optically thick model provides better 
consistence with the observations for the very broad component, while the optically thin 
scenario poses many problems. Hu et al. (2008) recently proposed another scenario that the 
classical BLR produces the H^vb component, meanwhile the H^b component is produced 
by an inflow of material. 

Comparing our new spectra with the previous ones allows us to identify a strong re- 
sponsivity to the continuum for the H/3vb component. As the continuum level at 5100A 
decreases from 5.2 x 10~^^ ergs cm^^ s^^ A^^ to ~ 2 x 10~^^ ergs cm~^ s~^ A^^, the mea- 
sured flux of H/3vB is reduced by a factor of 2-3. The corresponding equivalent widths (EWs) 
maintain a constant value ~ 100 A. All the EWs are measured against the continuum levels 
at 5 100 A. As an additional investigation, the upper panel in Figure 3 shows the smoothed 
fine profiles for Ha (the red lines) and H/9 (the blue fines). The profiles are produced by 
a subtraction of the modeled narrow-line profiles, and are smoothed by a box of 10 pixels. 
We present the average Ra/R/3 ratios as a function of radial velocity in the lower panel. 
The average values and the corresponding errors are calculated after binning the spectra 
over 1000km s~^. The results from the 2010 and 2011 spectrum are plotted by the red and 
blue symbols, respectively. In both observations, the calculated Balmer decrements have a 
plateau with Ha/H/? ~ 3 from -5000km s~^ to -|-5000km s~^, and decrease rapidly in the 
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wings. The models in Snedden & Gaskell (2007) indicate that the decrements at the wings 
can not be reproduced by the optically thin photoionization model in which the predicted 
Balmer decrement is rather uniformly around the value of Case B recombination. Snedden 
& Gaskell (2007) claimed that the decreased line ratios at wings in their sample are best 
explained by the optically thick model. 

The object shows peculiar variation for the H/3b component. When compared with the 
1990 spectrum, the component is decreased by a factor of 10 in flux in the 2000 spectrum. 
However, the flux returns to the 1990 level in the 2010/2011 spectra in which the continuum 
level is comparable with that of the 2000 spectrum. In addition to the flux change, the 
2010/2011 spectra show a large variation in the bulk blucshift velocity and line width for 
the EC when compared with the previous observations. The velocity and line width in the 
2010 (2011) spectrum are larger than those of the 1990 (2000) spectrum by nearly an order 
of magnitude. Putting the variations in flux and in velocity together suggests a possibility 
that the H/3b line might be emitted from an episodic outflow from the central engine. The 
outflow is at flrst radially accelerated at the launch point, which results in the increasing 
bulk blueshift and line width. The outflow then reaches a terminal velocity at larger distance 
where the radiative force and gravitation are negligible, and has a reduced optical depth. 
The dynamical timescale of BLR is Tdyn = -Rblr/A^, where i?BLR is the distance of the 
BLR from the central SMBH and AV ~ 10^ km s~^ is the typical velocity. Estimating 
the distance as -Rrlr = 20(L5ioo/10^^ erg s~^)°-^''' Id (Kaspi et al. 2000) yields a dynamical 
timescale Tdyn ~ lOyears, which does not pose fundamental challenge for the above scenario. 
The second possibility is that the object is a "Double-peaked" emission-line QSO whose 
broad Balmer line proflles are produced by an accretion disk (e.g., Chen et al. 1989; Gaskell 
2010). In contrast with the object, the typical "Double-peaked" emitters show line proflles 
with two peaks connected by a depressed plateau at the line center, and show rapid line 
proflle variation with a timescale typical of a few months. Finally, the current data can 
not exclude a possibility that the two Balmer components come from two separated SMBHs 
orbiting with their BLRs within a common NLR (e.g., SDSS J153636.22+044127.0, Boroson 
& Lauer 2009; Lauer & Boroson 2009). Unusual broad line profile (e.g., two displaced peaks 
or large shifted broad fine) is usually used as an indicator for binary SMBHs candidates (see 
recent review in Popovic 2011 and the references therein). The binary BLR model indicates 
a significant variation of line profile for different BLR orbiting phases. In PG 1416-129, 
assuming a Mbh ~ 10^ for each of the binary SMBHs, the Eq. (2) in Popovic (2011) 
results in a distance between the two black holes of ~ 10^^ pc when the orbital period 
of ~ lOyears is used. Alternatively, the signiflcant variability could be produced from the 
merging of the BLRs when the two black holes are in the "sim-detached" configuration. 
Long term spectroscopic monitor, especially in X-ray, is required to test which hypothesis is 
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correct for the object. 

4. CONCLUSION 

The spectra of luminous radio-quite quasar PG 1416-129 were re-observed in 2010 and 
2011 to study the variation of its broad-hne emission. The two spectra show neghgible hne 
variation within a timescale of one year. By re-analyzing the previous spectra, we identify a 
strong response to the continuum level for the very broad component. The wide wavelength 
coverage of the new spectra indicate a flat Balmer decrements at the line wings. With 
these results, we argue that a separate inner optically thin emission line region might not be 
necessary in the object. 
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Fig. 1. — The four solid black curves from top to bottom show the observed spectra of 
PG 1416-129 taken on 1990 February 17, 2000 May 08, 2011 May 29, and 2010 January 
25. For the latter two spectra, both A- and B-band telluric absorption features at A7600 
and A6800 are removed from the observed spectrum by the standards. The plotted spectra 
are vertical shifted by an arbitrary amount with respect to the 2010 spectrum for visibility. 
The shifted amount is 1.5 x 10~^^ erg s~^ cm~^ 1-2 x 10~^^ erg s~^ cm~^ and 

5.0 X 10~^^ erg s~^ cm~^ A"^ for the 1992, 2000, and 2011 spectrum, respectively. The 
modeled continuum and Fe II complex are overplotted by the red curve for each spectrum. 
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Fig. 2. — Emission-line profile modelings by a sum of a set of Gaussian components for 
the 1990 spectrum (top-left panel), the 2000 spectrum (top-right panel), the 2010 spectrum 
(middle panels), and the 2011 spectrum (bottom panels), after the continuum is subtracted 
from each of the observed spectra. For the latter two spectra, the left panels show the 
fitting for the H/3 region, and right panels for the Ha region. In each panel, the observed 
and modeled fine profiles are plotted by the fight and heavy solid fines, respectively. Each 
Gaussian component is shown by a dashed line. The sub-panel underneath each emission-line 
spectrum illustrates the residuals between the observed and modeled profiles. 
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Fig. 3. — Upper panel: Smoothed Ha (the red hnes) and H/3 (the blue hnes) hne profiles. 
The profiles from the 2010 spectrum are shown by the solid and dashed lines. The dotted and 
dot-dashed lines show the profiles from the 2011 spectrum. Lower panel: Average Balmer 
decrements across the line profiles for the 2010 (red symbols) and 2011 (blue symbols) 
spectra. The values of Ha/H/3 ratio and the corresponding errors are calculated by binning 
the spectra over 1000km s~^. 
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Ta ble 1: EMISSION-LINE MEASUREMENTS FOR PG 1416-129 

1990 FEBRUARY 17^ 



Flux FWHM AV 

LINE IDENTIFICATION IQ-^^ erg s"^ cm-^ km s"^ km s'^ 



H/3b 27.8 ± 1.4 4483 ± 194 -185 ± 50 

H/3vB 53.6 ±2.6 10963.6 ± 958 1297 ± 541 

H/3b+vb 81.4 ±3.0 

2000 MAY 08^ 



Flux FWHM AV 

LINE IDENTIFICATION 10-^^ erg s'^ cir-^ km s'^ km s"^ 

H^B 1.9 ±0.8 1443.4 ±441.3 308.0 ± 203. 

H/3vB 23.2 ±0.9 8548.1 ±432.7 1244.1 ±227.3 

H/3b+vb 25.1 ±1.2 

2010 JANUARY 25 ^ 



Flux FWHM AV 

LINE IDENTIFICATION 10"^^ erg s"^ cm-^ km s"^ km s'^ 

RpB 18.9 ± 1.2 4983 ± 248 -1356 ± 184 

H^VB 19.2 ± 1.0 10329 ±393 2393 ± 191 

H^B+vB 38.1 ±1.6 

HftB 63.8 ±4.6 4371 ± 166 -1160 ± 100 

HavB 94.2 ±4.1 9301 ± 189 2096 ± 160 

Hc^B+VB 158.0 ±6.2 

2011 MARCH 29 



Flux FWHM AV 

LINE IDENTIFICATION 10"^^ erg s"^ cm-^ km s"^ km s"^ 

H^B 16.6 ±2.2 4614 ± 356 -1059 ± 192 

H^VB 16.0 ±4.5 9117 ±2076 2684 ±486 

H/3b+vb 32.1 ±1.6 

HftB 59.8 ± 5.4 4397 ±210 -1642 ± 142 

HavB 106.0 ±9.4 7183 ± 163 1200 ± 150 

HaB+VB 165.8 ± 10.8 



"Specific flux - = 5.2 x lO^^^ ergs cm-^ A-\ 

''Specific flux - = 2.0 x lO^^*^ ergs cm-^ A'^. 

'^Specific flux F ° = 2.3 x 10"^*^ ergs cm^^ A'^. 

''Specific flux F - = 3.1 x IQ-^'^ gj.gs cm-^ s'^ A'^. 

^ 5100A ° 



